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Item topic Of Limite amplitude acoustic standing waves in 
sloshing modes of a cylindrical cavity is interesting 
maeeretically, but development of the subject has not been 
extensive. The purpose of this research is to present the 


meeults Of a power series perturbation approach to the 


peoo Lem . 





II. THE NON-LINEAR WAVE EQUATION 


A. GENERAL 

It is well known | that. for and kG nieve 
measures the fractional loss per wavelength and M is the peak 
Mach number of the source, loss terms and nonlinear terms in 
the constitutive equations can be separately approximated with 
the help of linear, lossless acoustic relations. The force 
equation appropriate for acoustical processes in systems for 


weich gravitational effects are unimportant is 


:; +(U-Y)U+ sup <eLU Oa 


a 
mee =ULlis Gees pateurcle vyeloer ty , Cis the instantaneous 
Sem@etty Of the fluid, P is the instantaneous total pressure in 
meet ield, and the operator J symbolically describes those 
physical processes leading to absorption and dispersion. We 


Mean two additional equations.:' The first is the equation of 


meome fOr a perfect gas 
P “Ses IK (Qeiaa) 


where r is a constant whose value depends on the particular 
gas involved and {k is the absolute temperature in Kelvin. 


igemeecond is the continuity equation 
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where s = ( F-F) / G is the condensation at any point and 4 
momene equilibrium density of the fluid. If we ignore 


rotational effects, then 


U =V¢ a 


where e is the velocity potential. Combination of Eqs.(2-1)- 
(2-2) and the neglect of terms of orders higher then Me, MA /f) , 


and (ALR) Z, yields a quadratically nonlinear wave equation, 
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WY is the ratio of heat capacities, and 
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(2-4) 


Mijemiert-hand side of Eq.(€2-3) is the classical, linear wave 


equation with losses pertinent to the system under study. 


The right-hand side can be interpreted as a forcing function 
consisting of a three-dimensional spatial distribution of phase- 
Mmemereit Sources. 

In a second-order perturbation theory, this volume forcing 
function is obtained from the classical(first-order) solution 
faeeere che acoustical problem. The second-order perturbation 
Peron Pa+P, describes the non-linearities P> resulting from 
Beem t interaction of the classical solution Pj. 

Peesiter—-Oorder perturbation solutions consider the interaction 
Secene nonlinear solution with itself, and the forcing function 
is composed of products of both classical and non-linearly 
Mew@eraced terms. Thus, if a system is driven at frequency w, 
Em@emmon—-linear term in Eq.(2-3) will force the existence of 

all integer multiples WW of the driving frequeney and the full 
solution must contain ali harmonics of the input frequency. 

in a closed cavity, each of those nonlinearly generated waves 
whose frequency lies near the resonance frequency of a standing 
wave of the cavity and whose assciated spatial function matches 
that of the standing wave can be strongly excited [ 2] : 

As far as the author has been able to determine, there has 
been only one previous study of this system published in the 
open literature. This was by Maslen and Moore in 1956 [3] 

Their approach resulted in a series expansion which did not 
converge if the relevant nomal mode frequencies were integerally 
related. Their interpretation of the quenching of the nonlinear 


effect was based on the "scattering effect of the wall". 





We feel our interpretation based upon nonlinearally generated 
volume sources stimulating the allowed standing waves of the 
cavity is more accessible and informative. Further, the 
mathematical approach developed herein appear to avoid any 
Sreeraeuities im convergence. Their conclusion that high 
amplitude monofrequency transverse oscillations can exist is 


@emsistent with our finding. 


Peer PLICATION TO THE CYLINDRICAL CAVITY 

Weewctrcular cylinder 1s provided with a "point" source of 
sound. By properly positioning the cylinder with respect to 
ieemoound Source if is possible to effectively drive the 
enclosed air into various modes of vibration. The rigid- 
walled cylinder forces the component of particle velocity 
perpendicular to each cavity surface to vanish at the surface. 
The resulting steady-state solution to the linear wave equation 


mmeeylindrical coordinates is 


P er Q, x)= Anmo’s Re I) CoS n§) COS (Wang) Jon Ram Y) 


(2-5) 


iemems are the cylindrical Bessel functions and application 


of the boundary condition to the sides yields 


4 — Jon (2-5-a) 
Yavin 


10 





were a is the radius of cylinder and Hem are tnewamsuments 
th | 

ef the exterema of the “‘J}/”" Bessel function. 

The nomal mode frequency is dependent on “Ray and ZA : 
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The standing wave well be identified by the ordered integers 
Cn. mm. 2U ) describing its spatial dependence. 
1. Symmetric Modes 
The simplest waves in cylindrical coordinates are those 
that depend only on the distance r from Z-axis, the gradient 


fees the form 
a ae 
af _ 4 a4 (2-6) 


where 


U= RY 


(2-6-a) 


and k=w/c is the wave number or propagation constant. 


The Laplacian becomes 
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The suitable trial solution appropriate for symmetric modes is 


the linear solution(0, m, 0) of frequency w 


Ia 








= = I. C4) AUL we (2-7) 


The velocity potential for this limit can be defined as 


= Mae (4) CoS WL (2-8) 
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and eC, 1s 


4 = 4 ™ J. C4) ces wl =-4MJ.4¢ wt (2-9) 


Because we assume that the surfaces of the cavity are rigid, 
S@emeae =a the appropriate boundary condition is J,'=0. 

vs 
Thus, we have ail Jian : 


To generate the second order solution, we first note that 





( v3) = = M Je (! = te AW ) (2-10) 
and 
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Meewen tne second derivatives of Eqs.(2-10) and (2-11) with 


respect to time are 
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Pi@eeertuting Eqs.(2-2) and (2-13) into the first term of the 


right-hand side of Eq.(2-3) yields 


zWy 
meee | ne AW | c. Ye *] (2-14) 
and the second term of right-hand side of £q.(2-3) becomes 
Zz a re 
WITl(-+4J =3),) (J.-J J. oes mot] C= 15) 
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x x xz Pa z x 
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re x = | z x —— Zz 
c. \V/ ar =ziw (> J. i ilk = ) (2-17) 


With the use of Eqs.(2-14) and (2-15) we can write the right- 


its 





Bend side of Eq.(2-3) as 


and the appropriate inhomogeneous wave equation for the second 


order perturbation solution is 
ee ; = 
rota yis — GW = +3 -2 J+(-21 
(CLI +H4 )ee = 4 My gle Ti WF bo JT 4 
xX 
+J-O-3) 1 Joes ZWT Qe) 


2. Non-Svmmetric Modes 
The non planar waves in cylindrical coordinates are 
those that depend on the distance r from Z axis and the angle 


§) from X-axis. The gradient takes the form 


a mee 


Y= (4549955) C=20)) 


q= kr 
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and the Laplacian becomes 
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Mremcuditabple trial solution appropriate to the forcing function 


of frequency w for exciting the (n, m, 0) standing wave is 
./ | 
al = Var i estas Cael WU ; jam / 0 (2-21) 
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where Tn (RY) is the Bessel function of the first kind and 


order n. The velocity potential is approximated by 


} = ee} | Tn CRY) CoS ng cos WT (2-22) 


Zig Ui/ co, 1s 


—— 


A , Ab 
Uae MCY Ty 25-79 cos wt —G “yf Sn Alt 1G cS WE ) G23) 
Co 


Thus, 


me 
(5) = MJ, wos TE AWN wl (2-24) 
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em@@mene second derivative of Eqs.(2-25) and (2-26) with respect 
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Smet a touting Eqs.(2-26) and (2-27) into the first term of 


Ea@m« 2-3) yields 
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In order to solve the second term of Eq.(€2-3), we can use 
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Mime oubstituting Eqs.(2-33) through(2-35) into the second 


memmeot EFq.(2-3) yields 
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With Eqs.(€2-30), (2+38) and (2-39) we can get the inhomogeneous 


wave equation for the second order perturbation as 
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meee Pert m=0 an Eqt2-40), it reduces to Eq.(2-19), as it must. 
BUbstituting n=1 into Eq.(2-40) yields the equation with 


forcing term resulting from the (1, m, 0) sloshing mode, 
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Recall that the equation with forcing term resulting from 
the(l1, m, 0) sloshing mode can be written as a function of 


PPecuency, 
aie = | ~ 
oe M cas xivt corse a oun 
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and the left-hand side of Ea.(2-3) is 
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If the harmonics of the frequency at which the cavity is driven 
peewmot Close to any resonant frequencies, we can ignore the 


Meesy term. So, Eq.(3-1-a4) can be written as 


(RG Hg)” 3) Re 7” 


Meus assume that the solution of Eq.(3-1) can be expressed as 
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and define 
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Pe@eenation of Eqs.(3-2) and (3-3) yields 


( 





Der st we, ay— lf | ++, _1+* 
sai) ap ty 3g) + J] Usp= af ib Sea Ge 


ae: a oe : oat ee 7" LTne| (3-6) 











= Ly) 7 
~ Ine E-2) S ai 





fee OWER SERIES METHOD 


In order to solve the right-hand side of Eq.(3-7), we can 
use the definition [4 
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Eq. (3-44) must equal the right-hand side of Eq.(3-40). 
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Recall that we can Beit UE ee Bde. Gs erneoucn(s= 39) 
eumemeeae. (3-14), (3-17), (3-20). Thus, with the use of Eq.(3-31), 
we can get V.. From Eqs.(3-46) through(3-48) and Eqs.(3-11), 
e720). (3-23), (3-26), (3-29), we can compute [ny . Pigs owl ea 
the use of Eq.(3-14), we can get U,. Therefore we can compute 
foneemend Ujty). Let us calculate the finite amplitude effects 
resulting from the nonlinear distortion of a forced radial mode. 
If we excite a(0, m, 0) mode and obtain the pressure at the 
circumference, then 4 = fogs Jo Given this value of +] , use 
emmpeags. (3-14), (3-17), (3-20) give the quantities C, D, E. 
From Eqs.(3-35) and (3-36) we can compute One . Thus, we can get 
Vasy) from Eq.(3-31). Let us calculate the finite amplitude 
effects resulting from the nonlinear distortion of a forced 
Sloshing mode. If we excite a(l, m, 0) mode and obtain the pre- 
ssure at the circumference, then 4 = HQ = hm . Given this value 
of Meise of Eqs.(3-11), (3-20), (3=23), (3-26), (3-29) give the 
quantities B. E, G, H, 0. From Eqs.(3-45) we can compute bade 
tities we Can get U,(y) from Eq.(3-41). Substituting n=0 at 
radial modes into V,(y) and U,(y) yields Vgly)=Ugj(y). Therefore, 
Reece use Of Eq.(3-3), we can get the solution of Eq.(3-2). 
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